Description of light clusters in relativistic nuclear models 
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Light clusters are included in the equation of state of nuclear matter within the relativistic mean 
field theory. The effect of the cluster-meson coupling constants on the dissolution density is dis- 
cussed. Theoretical and experimental constraints are used to fix the cluster-meson couplings. The 
relative light cluster fractions are calculated for asymmetric matter in chemical equilibrium at finite 
temperature. It is found that above T — 5 MeV deuterons and tritons are the clusters in larger 
abundances. The results do not depend strongly on the relativistic mean field interaction chosen. 
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I. INTRODUCTION 

To gain a deeper understanding of the physics involved 
in stellar core collapse, supernova explosion and pro- 
toneutron star evolution, an equation of state capable of 
describing matter ranging from very low densities to few 
times the saturation density and from zero temperature 
to a few MeV is needed. 

The crust of the star is essentially determined by the 
low density region of the equation of state (EOS), while 
the high density part will be important to define proper- 
ties such as the star mass and radius. Understanding the 
crust constitution of a neutron star is an important issue 
because it influences the cooling process of the star and 
plays a decisive role on quantities such as the neutrino 
cmissivity and gravitational waves emission. 

In the outer crust, where the density is lower, the for- 
mation of light clusters in nuclear matter will be energet- 
ically favorable at finite temperature, as in core-collapse 
supernovae or neutron star mergers, whereas in catalyzed 
cold neutron stars only heavy nuclei appear. At very low 
densities and moderate temperatures, the few body cor- 
relations are expected to become important and the sys- 
tem minimizes its free energy by forming light nuclei like 
deuterons (d = 2 H), tritons (t = 3 H), helions (h = 3 He) 
and a-particles ( 4 He) due to the increased entropy [l], H] . 
Eventually, these clusters will dissolve at higher densities 
due to Pauli blocking resulting in an homogeneous matter 
0. In particular, the cooling process of a protoneutron 
star is affected by the appearance of light clusters [|[ . 

The inclusion of light clusters (d, h, t and a parti- 
cles) in the nuclear matter EOS is discussed in [l|, where 
the most important thermodynamical quantities are cal- 
culated within a density-dependent relativistic model. 
The conditions for the liquid-gas phase transition are ob- 
tained, and it is seen how the binodal section is affected 
by the inclusion of these clusters. Moreover, an EOS is 
obtained starting at low densities with clusterized matter 
up to high density cluster-free homogeneous matter. In 
this work the density and temperature dependence of the 



in-medium binding energy of the clusters was determined 
within a quantum statistical approach Q and included 
in a phenomenological way in the relativistic mean field 
(RMF) model with density dependent couplings. 

The a particle is the most strongly bound system 
among all light nuclei and it certainly plays a role in nu- 
clear matter as has been pointed out in [HQ-Hi- Lattimer 
and Swesty [f| worked out the EOS in the compressible 
extended liquid drop model based on a non-relativistic 
framework appropriate for supernova simulations, for a 
wide range of densities, proton fractions and tempera- 
tures, including the contribution of a particle clusters. 
An excluded volume prescription is used to model the 
dissolution of a particles at high densities. The same is 
done by H. Shen et al. in Q where non-uniform matter 
composed of protons, neutrons, a particles and a single 
species of heavy nuclei is described with the Thomas- 
Fermi approximation and the TM1 parametrization of 
the non-linear Walecka model (NLWM). At low densi- 
ties, these particles arc described by classical gases. 

In 0] the virial expansion of low-density nuclear mat- 
ter composed by neutrons, protons, and a particles is 
presented, and it is shown that the predicted a particles 
concentrations differ from the predictions of the EOS pro- 
posed in [j| and [f| . The virial expansion was extended to 
include other light clusters with A < 4 @, Q . All possible 
light clusters were also included in recent EOS appropri- 
ate for simulations of core-collapse supernovae based in a 
nuclear statistical equilibrium model, including excluded 
volume effects for nuclei ficj . 

The aim of the present work is to study the disso- 
lution density of light clusters, i.e. the density above 
which light clusters do not exist anymore, in uniform 
nuclear matter within the framework of the NLWM [ll[ . 
The dependence of the dissolution density on the cluster- 
meson couplings as well as on the proton-fraction will be 
studied. In many models which include light clusters 
the dissolution density is determined by excluded vol- 
ume effects [5|, [(| [l(| . In Ref. @ , a statistical excluded 
volume model was compared with two quantum many- 
body models, a generalized relativistic mean-field model 
[l| and a quantum statistical model [U, [l2| ■ It was shown 
that the excluded volume description works reasonably 
well at high temperatures, partly due to a reduced Pauli 
blocking. However, at low temperatures the excluded 
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volume approach shows crucial deviations from the sta- 
tistical approach. 

In [13| . using a RMF approach, the coupling of the 
a-clusters to the a;-meson was employed to describe the 
dissolution. In this reference it was assumed that the in- 
tensity of the lu meson-cluster coupling was proportional 
to the mass number of the light cluster and no coupling 
to the cr-meson was included. 

The light clusters (d, h, t and a particles) are treated 
as point like particles (without internal structure) within 
the RMF approximation. Just as the nucleons in nuclear 
matter, the clusters interact through the exchange of er, 
uj and p meson fields. In order to understand how depen- 
dent are the results on the NLWM parametrization we 
consider three different parametrizations; NL3 [Til ], with 
a quite large symmetry energy and incompressibility at 
saturation and which was fitted in order to reproduce 
the ground state properties of both stable and unstable 
nuclei; FSU [H| and IU-FSU , which were accurately 
calibrated to simultaneously describe the GMR in 90 Zr 
and 208 Pb, and the IVGDR in 208 Pb and still reproduce 
ground-state obscrvables of stable and unstable nuclei. 
Furthermore, the IU-FSU was also refined to describe 
neutron stars with high masses (~ 2M Q ). 

In Ref . [l| , where a generalized RMF model with den- 
sity dependent couplings was developed including clus- 
ters as explicit degrees of freedom, the medium effects on 
the binding energy of the clusters were introduced explic- 
itly through density and temperature dependent terms. 
We aim at getting a description of the light clusters in 
nuclear matter that is simpler to implement, yet more 
realistic than the excluded volume mechanism. 

A brief review of the formalism is presented in section 
im the results and discussion of the meson-cluster cou- 
pling constants are given in section IIIII and conclusions 
are drawn in section ITVl 



II. THE FORMALISM 

We consider a system of protons and neutrons inter- 
acting with and through an isoscalar-scalar field tr with 
mass m Sl an isoscalar- vector field w M with mass m v , an 
isovector- vector field with mass m p and light clusters: 
d, h, t, and a particles. The Lagrangian density of this 
system reads: 

C ^ £>j ~t~ C a ~t~ &d ~\~ £u> £p ~\~ ^wpi 

j — p,n,i,/i 

where the Lagrangian Cj is 
with 



(1) 



iD 



$ = id" - g> vU V - ^ 



t ■ b^ , 



M* 



Mj fa, j=P,n,t,h 



(2) 
(3) 



where M = (M p + M n )/2 = 938.918695 Me V, M h = 
3M - B h , M t = 3M - B u M a = AM - B ai M d = 
2 M — Bd, with the binding energies given in Table [i] The 
a particles and the deuterons are described as in [l| 

£<* = \^K<t> a r(iD^ a )-^l(M* a ) 2 cp a , (4) 
£d = -(iD^-iD^)*(iD d ^ dv -iDfo^) 



where 



and 



(5) 

(G) 
(7) 



M* = M,-fa 
with j — a, d. The meson Lagrangian densities are 

Cup = Kglglu^hp, ■ b M 

where Q M „ = c^av-cU^, B M „ = d^-d^b^-gpib^ x 
h v ). The parameters of the models are: the masses and 
the coupling parameters <?<,, g v , g p of the mesons and A, 
k, £, and A„ of the nonlinear meson terms. In the above 
Lagrangian density r is the isospin operator. When the 
NL3 parametrization is used, £ and A„ are set equal to 
zero. The FSU and IU-FSU parametrizations include the 
nonlinear ujp term. 

The equations of motion in the relativistic mean field 
approximation for the meson fields are given by 



9 s Pi 



i=p,n,t,h 



ml^ + -A 9 t{^f+2K g l g y{blf = 



9vPi 
i—p,n,t,h,d,a 



m>° + 2A^(o; ) 2 &° = Yl S%Pi 

i—p^n,t : h 

where, p % s is the scalar density given by 



Ps = 



2S* + 1 ^,2,, Mi- fa 



2tt 2 



k z dk 



\Jh* + (Mi-faf 



Pi is the vector density, I\ (S l ) is the isospin (spin) of 
the specie i for i = n,p, t, h. At zero temperature all the 
a and deutcron populations will condense in a state of 
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zero momentum, therefore p a and pd correspond to the 
condensate density of each specie. 

The values of the parameters for each one of the rela- 
tivistic effective interactions, discussed below, are in Ta- 
ble lU and the corresponding properties of infinite nuclear 
matter at saturation density are shown in Table Hill 



TABLE I: Parameters for the cluster binding energy shifts 
defined in Eqs. ©- (JTOj) at T = MeV (taken from Q). 



The total energy density is given by 



Cluster i 


aj,i 




B9 [17J 




' MeV 5/2 fin 3 } 


(MeV) (MeV) 




t 


69516.2 


7.49232 


8.481798 


h 


58442.5 


6.07718 


7.718043 


a 


164371 


10.6701 


28.29566 


d 


38386.4 


22.5204 0.2223 


2.224566 



TABLE II: Parameter sets for the three models used in this 
work. 







NL3 [14] 


FSU [15] 


IU-FSU [16] 


m s 


(MeV) 


508.194 


491.500 


491.500 


id, 


(MeV) 


782.501 


782.500 


782.500 


m p 


(MeV) 


763.000 


763.000 


763.000 




9 2 s 


104.3871 


112.1996 


99.4266 




2 

9v 


165.5854 


205.5469 


169.8349 




9l 


79.6000 


138.4701 


184.6877 


K 


(MeV) 


3.8599 


1.4203 


3.3808 




A 


-0.015905 


+0.023762 


+0.000296 






0.00 


0.06 


0.03 






0.000 


0.030 


0.046 



TABLE III: Bulk parameters characterizing the behavior of 
infinite symmetric nuclear matter at saturation density: the 
saturation density po, the binding energy E/A, the incom- 
pressibility K, the symmetry energy a sym and its slope L, 
and the nucleon effective mass M* . 



NL3 FSU IU-FSU 



Po(ftn- J ) 
E/A{MeV) 
K ( MeV) 

(MeV) 
L ( MeV) 
M*/M 



0.148 

-16.24 

271.5 

37.29 

118.2 

0.60 



0.148 
-16.30 
230.0 
32.59 
60.5 
0.62 



0.155 
-16.40 
231.2 
31.30 
47.2 
0.62 




' kfdk^k? + (Mi - gfry 
+ ±-\g^ W - (-°) 4 - \ m l ( b °) 



3^3 



- Kglg 2 ^ (blf + {M>-9> + 9>°)K- 

i—a,d 

At zero temperature the free energy density T and the 
energy density coincide. 

We define the global proton fraction Y p as 
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P P 



where p = p p + p n +4: Pa + 2 P d + 3ph + 3pt is the total bary- 
onic density. Defining the mass fraction yi — Ai(pij p) of 
the various species we obtain 

112 1 

Y v = Vv + 2 Va + 2 Vd + 3 Vh + 3 Vt 

for the global proton fraction. 



A. Cluster binding energy 

In [l[ using a quantum statistical approach, the fol- 
lowing empirical quadratic form was proposed for the 
medium-dependent binding energy shift at zero temper- 
ature of the various clusters 



AB l (p,T = Q) 
where 



61340) + ^SB*(0) 



(8) 



A 

The total nucleon density is given by p = p P ot + and 

a>i,i 



SBi{T = 0) = 

for i = t,h, a, and 

SB d (T = 0) 



3/2 
h[2 



2a 2 a' 2 



(9) 



(10) 



for the deuteron. The parameters an, en 2 and 0^,3, 
taken from are listed for symmetrical nuclear mat- 
ter (Y p = 0.5) in Table H 

The density where a cluster becomes unbound, Pdis, 
referred as dissolution density, is given by 

Bi( Pdis ) = B° + ABi = 0. (11) 

The dissolution densities obtained for the light clusters 
considered at zero temperature and in symmetric nuclear 
matter are reported in the Table IIV1 
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TABLE IV: Densities at which the clusters become unbound 
given by equation (fTTI) for T = MeV. 



0.0035 



t 



I, 



d 



Pdis (fm~ 3 ) 0.00183175 0.00144835 0.00439226 0.00044815 



III. RESULTS AND DISCUSSIONS 

Our results will be reported in the present section. The 
calculations will all be performed at zero temperature ex- 
cept for Section [ill El where results at finite temperature 
will also be presented. 

In the following we will consider each type of light clus- 
ters separately in order to study the dissolution density 
of each of them in nuclear matter. 

We assume for the cluster couplings the form 



9- 



ijj9v 

al=Pj9s, j = t,h,d,a 
9p = 8i9 P , i = t,h, 

g\ = 0, i = d,a, 



(12) 



where the parameters T]j, f3j and Si are constants. Next 
we determine the dissolution density (pdis) as a function 
of these parameters and study the dependence of pdis 
on the global proton fraction (Y p ). Experimental and 
theoretical information will used to fix the parameters 
r)j, fa and 8i. 



A. Simple choice of the coupling parameters 

In order to have a feeling of the effect of including light 
clusters in the EOS, we will first fix the cluster coupling 
constants as g\ = 0, g l v = Aig v for i = t,h,d,a and 
9p = 9 P = 9 P by fixing 



Si 
Si 



ij , — 0, 
1, i = t, h. 



j = d, t, h, a, 



(13) 



0, 



d, a, 



in (11211 . This choice corresponds to the one used in 
[l~3l . l20| and may be justified in the following way: a) since 
we do not know how nuclear matter affects the binding 
energy of the light clusters we take the a coupling con- 
stants of the light clusters as zero (g\ = 0); b) we expect 
that the ui meson will play the role of an excluded vol- 
ume, i.e., at large densities will hinder the superposition 
of clusters by giving rise to a repulsive contribution. We 
take the coupling g\ = Aig v and consider that this would 
be of the order of magnitude expected; c) we consider 
that the coupling of the light clusters to the p meson is 
defined by the corresponding cluster isospin. 

The dissolution density at a given proton fraction and 
zero temperature is determined from the crossing be- 
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FIG. 1: (Color online) The free energy density minus the 
mass contribution at T = MeV as a function of the bary- 
onic density for Y p = 0.1 (a) and Y p = 0.5 (b) for the FSU 
parametrization, with choice (|13p for the cluster couplings. 
The EOS designed by the light cluster i has the largest pos- 
sible fraction of clusters i and if, necessary, also protons and 
neutrons. 



tween the proton-neutron free energy and the free en- 
ergy of matter with the largest number of clusters it is 
possible to form. This is the most stable matter with 
clusters. The stable phase is the one with the lowest free 
energy density. At low densities the equilibrium phase is 
the cluster phase and at large densities it is the proton- 
neutron matter. In fact, from the Fig. [T] it is seen that 
this transition is characterized by a negative curvature 
of the free energy density, and therefore, it does not cor- 
respond to a thermodynamically stable phase transition. 
The inclusion of larger clusters is necessary to realisti- 
cally describe this transition but this is not the aim of 
the present work. 

The free energy density minus the mass contribution, 
T = £ — Mp, is plotted in Fig. [T] as a function of the 
baryon density for different EOS, namely the proton- 
neutron EOS without clusters and the proton-neutron 
matter with the largest possible fraction of one type of 
light cluster for the proton fractions considered (Y p = 0.1 
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FIG. 2: (Color online) The dissolution density pd ia for T = MeV as a function of total proton fraction Y p for the various 
clusters, with choice (|13[1 for the cluster couplings: a (a), deuterons (b), helions (c) and tritons (d). In each graph the curves 
represent NL3, IU-FSU and FSU from top to bottom, respectively. 



and Y p = 0.5). The curves in Fig. [T]were obtained for 
FSU but for the other parametrizations the results are 
qualitatively the same. From Fig. [1] we see that there 
is some range of densities in which the free energy den- 
sity of the equations of state with one type of cluster is 
lower than the pure nucleonic matter (without any bound 
states) . In that range, the system can decrease its energy 
density by forming clusters. At T = 0, the EOS with a 
particles is by far the most stable followed by the EOS 
with tritons, helions and deuterons. This sequence of sta- 
bility is explained by the binding energy per nucleon of 
each cluster, having the a particle the highest value and 
the deuteron the lowest one. 

At finite temperature, all the clusters will be in chemi- 
cal equilibrium and their abundances will be determined 
by a balance between temperature, their mass and bind- 
ing energy. 

The calculated dissolution densities are shown in Fig. [2] 
for all the clusters and the three NLWM paramerizations 
with the choice (fT3")) for the cluster couplings. For the 
deuteron and a particles the maximum dissolution den- 
sity is at Y p = 0.5 allowing the conversion of all nuclcons 
into bound states (deuteron or a particles) . At Y p = 2/3 
and Y p = 1/3, all the nucleons are converted into he- 



lions and tritons, respectively, giving the highest dissolu- 
tion density. The NL3 and IU-FSU parametrizations give 
similar dissolution densities, with NL3 giving the highest 
ones. The IU-FSU paramctrization gives results between 
FSU and NL3, very close to NL3, therefore, it will not 
be considered in most of the discussion and only the NL3 
and FSU will be shown in the figures. The dissolution 
density of the light clusters within the different models 
is mainly determined by the isoscalar properties of the 
model. In particular, FSU has a quite large coupling 
constant to the a meson which defines a smaller effective 
nucleon mass at low densities and, therefore, clusters will 
dissolve at smaller densities. This is also the reason why 
the NL3 dissolution density is the largest: for the den- 
sities of interest NL3 has the smallest nucleon effective 
mass. 



We next generalize the choice (fT3"f varying smoothly 
one parameter (rji, Pi, Si) each time and keeping the 
other. 
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FIG. 3: (Color online) The dissolution density for tritons (a) 
and helions (b) as a function of 5 l (g p = 5i< 
(solid lines) and Y p = 0.1 (dashed lines). 



B. Varying the g p coupling constant 

To study the effect of the p meson-cluster coupling con- 
stants on the dissolution density for the hclion and 
triton, we fix the other coupling constants by g l v = Aig v 
and gl = and parametrize gl by g p = Sig p (i = h,t). 
The results are shown in Fig. [3l 

The densities of all species that are present are defined 
by the global proton fraction Y p . For values of Y p where 
only tritons and neutrons are present in nuclear matter, 
the free energy density term that depends on g p is 



9 P 
(m*) 



(S t pt + Pn) 



(14) 



and in the case where only helions and neutrons are 
present, 



1 9 ; 

^(Sh) = -- — {&hPh - PnY 



where (m*) = m 2 p + 2g 2 p g 2 A v (u; ) . 
Thus, for Y p = 0.1, an increase in g p through the pa- 
rameter Sh, reduces the symmetry energy of a system 
composed by helions and consequently, the dissolution 
density increases. For a system made of tritons the oppo- 
site happens, and an increase in g p increases the system 
symmetry energy leading to a decrease in the dissolution 
density. 

For Y p = 0.5, the helion and triton systems have a similar 
behavior, at some values of g l p and g p they acquire the 
lowest free energy density, and consequently, the highest 
dissolution density due to the vanishing of the symmetry 
term, HU) or ([TBI) . 

One important feature seen in Fig. |3]is the small range 
of variation of the dissolution density as a function of 
g p and g l p , compared with the range of variation of the 
dissolution density as a function of the io and a cluster 
coupling constants that we will discuss next. The cluster 
formation and dissolution is quite insensitive to the p 
meson-cluster coupling constant. 



C. Varying the gl coupling constant 

The effect of the w-cluster coupling constants on the 
dissolution density is studied by fixing the a and p cou- 
pling constants as gl — (i = t, h, d, a) and g p = g p 
(i = t,h). The gl is parametrized by rji as gl = r\ig v . 
The dissolution density versus r\i for all the clusters is 
shown in Fig. [3]for symmetric (Y p = 0.5) and asymmet- 
ric (Yp = 0.1) nuclear matter. 

An increase in g l v decreases the dissolution density for 
all clusters. This is due to the fact that an increase in gl 
leads to a more intense repulsion between the clusters, 
thus favoring their dissolution. 
For Y v 



(15) 



i p — 0.1 the dissolution density is lower than for 
t p — 0.5. This was expected because for Y p = 0.1 the 
number of clusters that can be formed is smaller than for 
symmetric nuclear matter. To dissolve a bigger number 
of clusters we need a bigger amount of energy and, thus, 
makes the cluster states with Y p = 0.5 more stable than 
Y„ = 0.1. 



D. Varying the gl coupling constant 

To study the dissolution density dependence on the 
cr-cluster coupling constants of all clusters g\ we fix the 
uj and p coupling constants as g p = g p (i = t, h) and 
gl = Aig v (i = t, h, d, a) and parametrize g\ as g\ = Pig s . 

The results are shown in Fig. [5] for Y p = 0.5 and 
ip = 0.1. As expected the dissolution density grows as 
the g l s increases. An increase in g l s makes the cluster 
more bound (smaller effective mass) and, therefore, clus- 
ter states become more stable. Except for the hclion, 
in all other cases for a given value of (3, the dissolution 
densities of the clusters become larger for Y p — 0.1 than 
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Y p = 0.5. This is because when considering Y p = 0.5 we 
are comparing matter with clusters with nuclear sym- 
metric matter, that is nuclear matter in its most bound 
state. For Y p = 0.1 nuclear matter is not even bounded. 
Increasing the binding energy of the clusters will stabi- 
lize them more relatively to the homogeneous asymmet- 
ric nuclear matter than in the previous case to symmetric 
matter. For hclions there are too few clusters in matter 
with Y p = 0.1 and therefore the curves Y p = 0.5 and 0.1 
will never cross. 

When fa is well below Ai, the dissolution density 
slightly changes for all clusters. For some value of 
below Ai the dissolution density increases abruptly and 
the system becomes absolutely stable, in other words, the 
clusters will not dissolve for any density. 



E. A physical choice for the cluster coupling 
constants 

After analyzing the effect of the different cluster-meson 
couplings, we discuss how these couplings could be fixed 
taking into account theoretical many-body calculations 
and experimental information. 

To fix the a and uj clusters couplings two constraints 



are required for each cluster. We consider the dissolu- 
tion densities for symmetric nuclear matter at T = 
MeV given in Ref. Q and shown in Table IIVI How- 
ever, this information is not enough to completely fix the 
cluster-meson couplings and we can find several sets of 
coupling constants that reproduce the dissolution densi- 
ties. Therefore, we choose to fix two coupling constants 
for each cluster and fit the remaining to reproduce its 
dissolution density, namely, we fix g\ and g l p . For g\ we 
take several values, 0, Ai/2g s , 3Ai/4g s , AAi/5g s , for all 
clusters and for g p we consider g l = g p , Si = 1, for triton 
and helion, and g l p = 0, <5; = 0, for the isospin symmet- 
ric clusters, which is a reasonable choice, and, as we have 
shown, this coupling does not affect much the dissolution 
density of the cluster. Then, for each choice of gl, gl is fit 
to reproduce the dissolution density of each cluster. The 
calculated values are shown in Table [V] Different choices 
of g % s , together with the above restriction for g' v , will al- 
low us to discuss the effect of this coupling constant on 
the overall particle fraction. In general we should have a 
different ratio gl/g s for each cluster, but, to simplify the 
discussion, we take the same for all clusters. 

The relative magnitude of the w— cluster coupling of 
the different models is related to the uj— nucleon coupling 
constant g v : the smaller the values of g v the larger the 
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FIG. 5: (Color online) Dissolution density of each cluster as a function of Pi (gl = j3ig s ) for Y p — 0.5 (solid lines) and Y p — 0.1 
(dashed lines): tritons (a), deuterons (b), helions (c) and a (d). 



to— cluster coupling fractions gl/g v , and the difference 
between gl for the models used is around 1%. In fact, 
for example, we find that for gl = 0, g\ is 43.22 (NL3), 
42.40 (FSU) and 42.92 (IU-FSU). 

More experimental information is needed in order to 
determine the coupling constants of the phenomenologi- 
cal model we propose. Recently, in-medium binding en- 
ergies and Mott points of d, t, h and a clusters have been 
determined experimentally [19j. For each cluster a com- 
bination of density and temperature corresponding to a 
vanishing in-medium cluster binding energy, the Mott 
point, were determined in a temperature range close to 
T = 5 McV. These results, together with the dissolution 
densities at T = MeV given in Table IIV1 will be used 
to determine the cluster coupling constant g\. 

In the study by Typel et al. |l| the binding energy of 
the clusters in the medium, 

Bi = AiM* - M* (16) 

where M* is the effective mass of cluster i defined in ^ 
and (J7J), was parametrized as a function of density and 
temperature. In Fig. |6ja) we compare the in-medium 
binding energy of the a-clusters as a function of the den- 
sity for several values of the couplings g\. All results 
were obtained at T = 5 MeV but, in fact, they do not 



depend much on the temperature. In the same figure we 
also include results obtained in [l| at T = 5 MeV, which 
turn out to be very close to the calculated binding en- 
ergies settings g\ = 3 A;/4 g s . In order to reproduce the 
results of [l| we must consider temperature dependent 
cr-cluster. We note, however, that choosing a tempera- 
ture dependent parameter will alter the usual thermody- 
namic relations between energy density and pressure. In 
Fig. [HIb) we compare the in-medium binding energy of 
the four light clusters obtained in our calculation setting 
gl = 3Ai/4 g s with the results of the best accordance 
occurs for the a-clusters. We may get a better accordance 
for all the clusters if a different coupling fraction gl/g s is 
taken for each cluster. 

In Fig. [7] we compare the Mott densities, the densi- 
ties at which the binding energies defined in (|16p van- 
ish, at T = 5 setting g* = 3Ai/4g s (open triangles') 
with the experimental Mott points determined in [19( 
(full dots with errorbars) and the parametrization of the 
Mott point given in |jj (dashed lines). Results in [l9[ 
are obtained at different temperatures for the four light 
clusters, however, all temperatures lie close to 5 MeV. 
We conclude that choosing g l s = 3Ai/4g s , and the cor- 
responding gl and g l p couplings defined in Table [Vj at 
T = 5 MeV in our model, gives results that are in rea- 



9 



1 Typel '2010' — 
g«/g s = 4A/5—- 

3 A /4 

A/2 -■- 






0.002 0.004 0.006 0.008 
p(fm- 3 



i • 1 • r 

Typel 2010 

gi/g s = 3Ai/i---- 




0.002 0.004 0.006 0.008 
p (fm" 3 ) 



TABLE V: Cluster coupling constants gl obtained for sym- 
metric matter (Y p — 0.5) that reproduce the dissolution 
densities of Table IIVI for several values of the coupling g\ . 
9p/9p ~ 1 f° r the triton and the helion and zero for the 
deuteron and the a. 

gl/gv 



gl/gs 





At/2 


2Ai/3 


SAi/4 


iAi/5 


FSU 












t 


2.961 


3.793 


4.177 


4.382 


4.509 


h 


3.283 


4.064 


4.428 


4.624 


4.745 


d 


3.038 


3.260 


3.422 


3.516 


3.577 


o 


4.440 


5.027 


5.439 


5.675 


5.824 


NT ^ 












t 


3.357 


4.172 


4.559 


4.765 


4.892 


h 


3.718 


4.485 


4.847 


5.043 


5.164 


d 


3.371 


3.585 


3.744 


3.837 


3.897 





4.919 


5.491 


5.897 


6.131 


6.278 


IUFSU 












t 


3.303 


4.126 


4.508 


4.715 


4.842 


h 


3.658 


4.429 


4.792 


4.988 


5.110 


d 


3.327 


3.543 


3.702 


3.796 


3.859 


a 


4.854 


5.429 


5.837 


6.071 


6.220 



FIG. 6: (Color online) Binding energy, as defined in Eq. (|16l) . 
of light clusters in equilibrium with symmetric nuclear mat- 
ter described within FSU at T — 5 MeV. In (a) the binding 
energy of Q-particles is given for several values of the cou- 
pling g sa (dashed, dotted and dash-dotted curves). The a in- 
medium binding energy given in Typel 2010 Q for T = 5 MeV 
is also included (full line). In (b) we consider gl = 3Ai/4g a 
for all the clusters (dashed lines) and compare with the corre- 
sponding results of Typel 2010 [J at T = 5 MeV (full lines). 
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FIG. 7: (Color online) The Mott density for different temper- 
atures and light clusters: results from Typel 2010 [l[ (dashed 
lines), experimental prediction by Hagel 2012 [li| (full dots 
with errorbars), results for gl/g s = 3A;/4 and T = 5 MeV 
(open triangle). 



sonablc agreement with the data of Ref. [T9| for the Mott 
density. 

As a test of the parametrizations proposed for g\ = 
3 Ai /Ag s in Table [V] we calculate the particle fraction 
at finite temperature and in chemical equilibrium, ^ = 
Zifip + (At - Zi)[i n with i = a,h,t,d pl|20|. We first 
focus on the effect of g l s and show in Fig. |5{a) the frac- 
tion of light clusters at T = 5 for symmetric nuclear 
matter obtained for two values of the er-cluster coupling 
constant, (dashed lines) and 3Ai/4g s (full lines). It 
is seen that the particle fraction and dissolution density 
is sensitive to this coupling: for a larger coupling, both 
the fraction of particles and the dissolution density are 
larger. In Figs. Mp) — Efci), we always set gl = 3Ai/Ag s . 
In Fig. [8jb) we compare two models, NL3 (dashed lines) 
and FSU (full lines) . It is seen that the results do not dif- 
fer much except for the deuterons that dissolve at a larger 
density and attain a larger fraction with NL3. At very 
low densities when the interaction between nucleons and 
clusters is negligible both models coincide; at densities 
just below the dissolution density, where models differ 
the most, the differences are not very large. The effect 
of temperature and isopsin asymmetry is shown in Fig. 
|S{c) and (d). In Fig. [Sfc) particle fractions in symmetric 
matter are compared for T — 10 MeV (full lines) and 
T = 5 MeV (dashed lines). At T = 10 MeV the deuteron 
fraction is already the largest fraction and the a fraction 
the smallest. In neutron rich matter, such as the one rep- 
resented in[5Ud), the deuteron fraction is the largest and 
the tritons come in second both at T = 5 and 10 MeV. 
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FIG. 8: (Color online) Fraction of light clusters in equilibrium with nuclear matter: (a) for FSU with Y p — 0.5 and T = 5 
MeV, and gl = (dashed lines) and gl = 3Ai/4g a (full lines); (b) for gl = 3Ai/4g a with Y p = 0.5 and T = 5 MeV and FSU 
(full lines) and NL3 (dashed lines); (c) for FSU with gl = 3Ai/4g s and Y p = 0.5, and T = 10 MeV (full lines) and T — 5 MeV 
(dashed lines); (d) for FSU with gl = 3 Ai/Ag s and Y p = 0.3, and T = 10 MeV (full lines) and T = 5 MeV (dashed lines). 



The lar ge a mount of deuterons is in agreement with re- 
sults of |18| , where composition of low-density supernova 
matter composed of light nuclei was calculated within a 
quasiparticle gas model, and 3, where composition of 
the outer layers of a protoneutron star was studied. 



IV. CONCLUSIONS 

In the present work we have included light clusters 
in the EOS of nuclear matter within the framework of 
rclativistic mean field models with constant coupling 
constants. Clusters are considered point-like particles 
that interact with the nucleons through the exchange of 
mesons. The formation of clusters is favorable at a quite 
low density; below 0.002 fin -3 , we expect that it is a 
reasonable approximation to consider them as point like. 

It was shown that the dissolution of clusters is mainly 
determined by the isoscalar part of the EOS. Therefore, 
the p-cluster coupling constant was simply determined 
by the isospin of the cluster. To fix the a and ui clusters 
couplings two constraints are required for each cluster. 
In the present work we have considered the dissolution 



density obtained in a quantum statistical approach Q as 
a constraint to fix the w-cluster coupling. Recent experi- 
mental results for the in-medium binding energy of light 
clusters [l9[ allow the determination of the er-cluster cou- 
pling strength at T ~ 5 MeV. The virial expansion of the 
EOS at low densities and finite temperature is another 
constraint that could be used to fix the cr-meson coupling 
and that will be investigated. 

We have applied the couplings proposed within the 
present to study symmetric and asymmetric nuclear mat- 
ter with light clusters in chemical equilibrium at finite 
temperature and we have determined the relative light 
cluster fractions at T = 5 MeV and T = 10 MeV. It was 
shown that a larger er-clustcr coupling gives rise to larger 
dissolution densities and larger particle fractions. 

The experimental determination of Mott points at 
more temperatures than the ones obtained in [19j would 
allow the determination of the temperature dependence. 
A comparison of the in-medium binding energies ob- 
tained within the present model with the ones proposed 
in [l| , indicate that the last may be reproduced if temper- 
ature dependent meson-cluster couplings are obtained. 
It was shown that deuterons and tritons are the clus- 



11 



ters with larger abundances in asymmetric matter above 
T = 5 MeV.. The results do not dependent much on the 
RMF interaction chosen. 
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